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Abstract — Energy harvesting has emerged as a powerful tech- 
nology for complementing current battery-powered communi- 
cation systems in order to extend their lifetime. In this paper 
a general framework is introduced for the optimization of 
communication systems in which the transmitter is able to 
harvest energy from its environment. Assuming that the energy 
arrival process is known non-causally at the transmitter, the 
structure of the optimal transmission scheme, which maximizes 
the amount of transmitted data by a given deadline, is identified. 
Our framework includes models with continuous energy arrival 
as well as battery constraints. A battery that suffers from energy 
leakage is studied further, and the optimal transmission scheme 
is characterized for a constant leakage rate. 

Index Terms — Battery leakage, battery size constraint, broad- 
cast channel, continuous energy arrival, energy efficient com- 
munications, energy harvesting, rechargeable wireless networks, 
throughput maximization. 

I. Introduction 

Energy efficiency is a key challenge in the sustainable 
deployment of battery-powered communication systems. Ap- 
plications such as wireless sensor networks depend critically 
on the lifetime of individual sensors, whose batteries are 
limited due to physical constraints as well as cost consid- 
erations. Power management is essential in optimizing the 
energy efficiency of these systems in order to get the most 
out of the available limited energy in the battery. A com- 
plementary approach has recently been made possible by 
introducing rechargeable batteries that can harvest energy from 
the environment. Several different technologies have been 
proposed and implemented for harvesting ambient energy such 
as solar, radio-frequency, thermoelectric or solar (see (Ti- 
lls) and references therein for various examples of energy 
harvesting technology). 

Harvesting energy from the environment is an important 
alternative to battery-run devices to extend their lifetime. How- 
ever, it is important to design the system operation based on 
the energy harvesting process to increase the efficiency. Energy 
harvesting systems have received a lot of recent attention ||6l- 
iflOl . Node and system level optimization have been considered 
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from both practical and theoretical perspectives. The previous 
work that are most relevant to the problems studied in this 
paper are |[TTI - lfT4ll . In ifTTI . the problem of transmission time 
minimization is studied when the data and the energy arrives 
at the transmitter in packets; and the transmission power is 
optimized when the data and energy arrival times and amounts 
are known in advance. In fl2l . the amount of transmitted data 
is maximized for an energy harvesting system under deadline 
and finite battery capacity constraints. Reference [12| also 
shows that the transmission time minimization problem studied 
in IfTTI and the transmitted data maximization problem are 
duals of each other and their solutions are identical for the 
same parameters. The problem is extended to the broadcast 
channel in |[T3l - lfT6l . to the relay channel in 1(171 . and to the 
multiple access channel in lfl8ll . 

The problem considered in this work is that of maximizing 
the amount of data that is transmitted within a given deadline 
constraint under various assumptions regarding the energy 
harvesting model as well as the battery limitations. Our 
focus is on the offline optimization of the energy harvesting 
communication system, that is, we assume that the energy 
arrival profile is known in advance. We first introduce a general 
framework for transmitted data maximization by adjusting the 
transmit power in an energy harvesting system with battery 
limitations. Our model includes continuous energy harvesting, 
generalizing the packetized energy arrival model considered in 
ifTTIl and lfl2l . Moreover, different from the previous work, our 
model also includes the realistic scenario of battery degrada- 
tion over time by considering a time-varying battery capacity. 
We show that these constraints can be modeled through 
cumulative harvested energy and minimum energy curves, 
which are then used to obtain the optimal transmission policy. 
The framework introduced for the energy harvesting system 
optimization is similar to the calculus approach introduced by 
Zafer and Modiano for energy-efficient data transmission in 
Ifl9ll . We later show that the proposed framework also applies 
to a broadcast channel with an energy harvesting transmitter. 

We then consider a more realistic battery model with energy 
leakage. Assuming a constant leakage rate, we identify the 
optimal transmission strategy for the case of a packetized 
energy arrival model. 

The paper is organized as follows. Section [Tj] presents 
the system model. Optimal transmission scheme for a point- 
to-point system under battery size constraints is derived in 
Sections [III] In Section IIV1 it is shown that the proposed 



framework can be used to characterize the optimal transmis- 
sion scheme in an energy-harvesting broadcast channel. We 
consider battery leakage in Section [V] and find the optimal 
transmission scheme for a linear leakage rate. Finally, conclu- 
sions are provided in Section [VI] 

II. System Model 

We consider a continuous-time model for both the harvested 
and the transmitted energy, that is, the harvested energy is 
modeled as a continuous-time process, while the transmitter is 
assumed to be able to adjust its transmission power, and hence, 
the transmission rate, instantaneously. This continuous-time 
model generalizes the discrete-time arrival model considered 
in ifTTI and lfl2l . A cumulative curve approach is used to 
described the flow of energy in the system. 

Definition 2.1 (Harvested Energy Curve): The harvested 
energy curve H(t) is a right-continuous function of time t, 
t, 6 M + , that denotes the amount of energy that has been 
harvested in the interval [0,i]. 

Definition 2.2 (Transmitted Energy Curve): The transmit- 
ted energy curve E(t) is a continuous function with bounded 
right derivative, that denotes the amount of energy that has 
been used for data transmission in the interval [0, t], t S K + . 

Naturally, we require E(t) < H(t), i.e., the transmitter 
cannot use more energy than that has arrived. We also consider 
a "minimum energy curve" that might model, for example, a 
battery size constraint. 

Definition 2.3 (Minimum Energy Curve): Given an har- 
vested energy curve H(t), a minimum energy curve M(t) 
is a function satisfying M(t) < H(t), Vt > 0, and denotes 
the minimum amount of energy that needs to be used by the 
system until time t. 

Given the harvested energy curve and the minimum energy 
curve, a feasible transmitted energy curve should satisfy the 
conditions M(t) < E(t) < H(t), Mt > 0. Among all feasible 
transmitted energy curves, our goal is to characterize the one 
that transmits the highest amount of data over a given finite 
time interval [0,T]. We consider offline optimization, that is, 
the harvested and the minimum energy curves are assumed to 
be known in advanced. 

We assume that the instantaneous transmission rate relates 
to the power of transmission at time t through a rate func- 
tion r( ), which is a non-negative strictly concave increasing 
function of the power with r(0) = 0. We note here that 
many common transmission models, such as the capacity of an 
additive white Gaussian noise channel, satisfy these conditions 
[ 19 1. The total transmitted data corresponding to a given curve 
E(i) over the interval [0, T] is found by 

V(E(t)) 4 [ T r (E'(t))dt, (1) 
Jo 

1 This is an accurate assumption for systems in which the energy harvesting 
process can be modeled as a deterministic process. For example, in solar based 
systems the amount of energy that can be harvested at various times of the 
day can be modeled quite accurately. In some other systems, harvested energy 
depends on the operating schedule of the harvesting device rather than the 
energy source, such as shoe-mounted piezoelectric devices; and the harvested 
energy curve can be modeled accurately in advance. 



where E'(t) is the derivative of function E(t) at time t, and it 
gives the power of transmission at that instant while r(E'(t)) 
is the corresponding transmission rate. 

III. Optimal Transmission Scheme under Battery 
Size Constraints 

In our problem formulation we assume that the transmitter 
always has data to transmit. Hence, the minimum energy 
curve can be used to model a constraint on the battery size, 
forcing the system to use any energy that cannot be stored 
in the battery for transmission of additional data before it is 
discarded. For a fixed energy curve E(t) and unlimited battery 
size, the energy that is available in the battery at time instant t 
is given by H(t) — E(t). However, if the battery size is b, we 
should have H(t) — E(t) < b. Consequently, the associated 
minimum energy curve is given by M (i) = ma,x{H(t) — b, 0}. 

We can also consider a time-varying battery capacity b{t), 
which can model the degradation in the battery capacity over 
time. This is a common phenomenon in rechargeable batteries 
used for energy harvesting applications. See Fig. |l(a)| for an 
illustration of the harvested and minimum energy curves for 
a battery with continuously decreasing capacity. 

Now, the optimization problem can be stated as follows. 

max T)(E(t)) = /f r(E'(t))dt (2) 

such that M(t) < E(t) < H(t),Vt G [0, T], (3) 

where T specifies the set of all non-decreasing, continuous 
functions with bounded right derivatives for all t G [0,T] and 
with E(0) = 0. 

We first present the optimality conditions for the transmitted 
energy curve. Similar to previous studies, such as ATI . fi"2l . 
lfl9l and [20], our main tool is the Jensen's inequality given 
in the following lemma (in the integral form). 

Lemma 3.1: [Jensen's inequality] Let / : [a, b] — > K, be 
a non-negative real valued function, and (/>(■) be a concave 
function on the real line, then 

with strict inequality if </>(•) is strictly concave, a ^ b, and / 
is not constant over the interval [a, b]. 

Consider the simple setup in which the battery has available 
energy Eo at time t = 0, no further energy is harvested, and the 
minimum energy curve is given as M{t) = for < t < T 
and M(T) = Eq. We will prove for this simple setting that 
the constant power curve transmits the maximum amount of 
data over the time interval [0,T]. 

For any transmitted energy curve E(t) with non-constant 
power, by replacing the function / in Lemma 13.11 with 
E'(t)/T, and letting a = 0, b = T and <f>(-) = r(-), we 
obtain 

{ f T E'(t) , \ f T r(E'(t)) , 
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(a) System with packet arrivals and a time-varying battery constraint. 
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(b) The "dying" battery example. 



Fig. 1. Illustration of the harvested and minimum energy curves for different examples. 



which is equivalent to 

(■T 



J r (E'{t))dt<Tr(^-Y (6) 

Note that Tr(y) is the transmitted data by the constant 
power scheme. Hence, this proves the fact that the maximum 
data is transmitted by this scheme. We can express this result 
in a more general context as in the following theorem. 

Theorem 3.2: Let E(t) be a feasible transmitted energy 
curve and S(t) be a straight line segment over interval [a, b] 
that joins E(a) and E(b), < a < b < T. If S{t) satisfies 
M(t) < S(t) < H(t) for a < t < b, the transmitted energy 
curve defined as 



E(t), te[0,a) 
E(t)={s(t), t€[a,b) 

E{t), te[b,T] 



(7) 



satisfies V{E(t)) > V(E(t)). 

The following theorems state, respectively, the uniqueness 
of the optimal transmitted energy curve and the optimality 
conditions. Their proofs follow similarly to those of Theorem 2 
and Lemmas 2-4 in |fl9l . 

Theorem 3.3: For a strictly concave rate function r(-), 
if E(t) is a feasible transmitted energy curve which does 
not have any two points that can be joined by a distinct 
feasible straight line, then E{t) is unique and it maximizes 
the transmitted data. 

Theorem 3.4: Let E opt (t) be the optimal energy expendi- 
ture curve and to be any point at which the power of transmis- 
sion changes, i.e., the slope of E opt (t) changes. Then, at to, 
E op t{t) intersects either H(t) or M(t). If E opt (t Q ) = H(t ), 
then the slope change must be positive. If E opt (to) = M(to), 
then the slope change must be negative. 

The optimal transmitted energy curve is also the one that 
has the minimum length, and hence, the same "string visu- 
alization" suggested in |19| can be applied here. The string 
visualization suggests that, if we tie one end of a string to 
the origin and connect it to the point (T,H(T)) tightly while 
constraining it to lie between H(t) and M(t), this string gives 
us the optimal energy expenditure policy. 
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Fig. 2. Continuous energy harvesting curve for a solar panel. 



A special case of the framework considered here is the 
one with packetized energy arrivals and without any battery 
constraint. This is the energy-harvesting dual of the packet 
arrival problem considered in [20|. As it is shown in 1 12 1, this 
is equivalent to the problem of transmission time minimization 
problem studied in 1111 . In this problem we have M(t) = 
for t G [0,T], and N energy packets arrive at times {ti}fL 1 . 
The algorithm that gives the optimal transmitted energy curve 
for this problem can be obtained following [11| and |fl9l . 

Another example that fits into the general structure intro- 
duced above is the following. Consider a wireless system with 
an energy storage unit consisting of N batteries. Assume that 
all the batteries are full initially and a total of £jv = ^TJili &i 
energy is available in the system at time t ~ 0, where bi 
is the capacity of battery i. It is assumed that the batteries 
in the system have finite lifetime, and they die at certain time 
instants, ti, i = 1, . . . , N. The problem is to find the maximum 
amount of data that can be transmitted until the last battery 
dies, i.e., until tjy. In this problem we have H(t) = Em 
for t G [0,£jv], and M{t) can be obtained as in Fig. |l(b)| 
Note that, since once the battery dies, the energy stored in 
it is not available for transmission anymore, and since we 
always have data in the queue to be transmitted, it is always 
beneficial to use the available energy in a battery before it 



dies. In this sense, we can consider the time until a battery 
dies as a deadline constraint on the time the available energy 
in this battery should be used. The optimal transmitted energy 
curve can then be found using the string argument as seen in 
Fig. [Kb)] 

As an example of continuous energy arrival, we consider 
here a model of a solar panel harvesting energy during the 
day. The amount of energy harvested per unit of time changes 
during the day. While no energy is harvested when there is no 
sun, the harvested energy is maximized at noon (see J2T|). 
We model the rate of harvested energy with the function 
h(t) = 5 - ^(t - 12) 2 for 6 < t < 18, and h(t) = 
elsewhere, where t € [0, 24] denotes the time of the day 
(hours), such that H (t) = J* h{r)dT. The unit of energy 
depends on the solar panel characteristics. The corresponding 
harvested energy curve is depicted in Fig. [2] 

Assume that we want to maximize the amount of data that 
can be transmitted up to time t — 18, i.e., until the panel 
stops harvesting energy. Based on the above arguments, the 
optimal transmitted energy curve is identified as follows. First 
we draw a tangent to the harvested energy curve from the point 
(18, H (18)), and denote its intersection with the curve by A. 
The transmitted energy curve follows the harvested energy 
curve from the origin up to A 0. Afterwards, it follows the 
straight tangent line, i.e., it uses constant power transmission. 
Note that, while it is easy to prove the optimality of this 
strategy using Theorem 13. 31 the discrete energy arrival models 
studied in ifTTI and lfl2l do not apply here. 

IV. Optimal Broadcast Scheme with Battery 
Constraint 

In this section, we show that the general approach intro- 
duced in Section [III] can be instrumental in identifying the 
optimal transmission policy in a broadcast channel (BC) with 
an energy harvesting transmitter [ 13 1, [ 14 1. Consider the same 
energy harvesting model at the transmitter as before; however, 
now there are two receivers in the system, and the transmitter 
has independent data for each receiver. 

The BC problem is studied in [13] and [14J; however, 
the solutions in these papers are elaborated from the basics 
rederiving the behavior of the optimal transmission policy in 
the BC scenario. Here, we show that the general approach 
introduced in previous section for the point-to-point setting 
can be directly applied to the BC scenario as well. This 
approach allows to generalize the results in iTPTl and (14 1 to 
continuous energy arrivals, and introduce battery constraints 
in the problem formulation flT31 . fl6l . 

We consider an additive white Gaussian BC in which the 
signal received at receiver i is given by 



Y, — X + Z, 



1,2, 



(8) 



2 In practice, a continuous adaptation of the transmission rate is unrealistic 
due to the block structure of channel coding, and the finite number of 
modulation and coding modes available. However, such practical constraints 
are out of the scope of this paper. 



where X is the channel input of transmitter and Zi is the zero- 
mean Gaussian noise component with variance Ni. Without 
loss of generality, we assume that N 2 > N% > 0. Let 
Bi (t) denote the total number of bits transmitted to receiver 
i up to time t. Our goal is to maximize the weighted sum 
of transmitted bits by time T, niBi(T) + fi 2 B 2 (T) for some 

Ml, ^2 > 0. 

In the broadcast channel setting, the transmitter not only 
needs to identify the transmitted energy curve E(t), but also 
has to decide how to allocate the power among the two 
receivers at each time instant. Accordingly, we denote by p\ (t) 
and p 2 (t) the power allocated to each receiver at time t. The 
optimization problem can be written as follows. 



max 

Pl(t),P2(t)>0 



Mi J T n (t)dt + n 2 J T r 2 (t)dt (9) 



such that M(t) < J *Pi(r) + p 2 (r)dT < H(t),t e [0,T] 

We assume that the rate-power functions are operating on the 
boundary of the capacity region of the Gaussian BC: 



ri(i) = i log2 ( 1 + £|> 



(10) 



(ii) 



The considered optimization can be decoupled into two max- 
imization problems as follows: 



max 

s(t)er J0 

M(t)<E(t)<H{t) 



max fJ.ir 1 (t) + /i2r 2 (t) 

Pl(i),P2(«)>0 

b>i(t)+P2(*)=pM 



dt. 



(12) 

where we define p(t) = E'(t). 

First, we consider the maximization problem in between 
brackets in d!21 l. Defining /1 = j^, we can make the following 
observations on its solution^: 

• If fj, > ^p, no power is allocated to the first receiver, i.e. 
Pi = 0, independent of the total power. 

• If // < 1, no power is allocated to the second receiver, 
i.e. p 2 — 0, independent of the total power. 

• When 1 < fi < jZ, the optimal power allocation behaves 
as follows. If the available total power is below pth = 
N2 ~^™ 1 , all the total power is allocated to receiver 1, 
i.e., pi — p and p 2 = 0. On the other hand, if p > p t h, 
then we have p x = p th and p 2 = p - pth- 

Note that, if /j, > ^ or fi < 1, the problem reduces to the 
point-to-point setting; hence, we assume 1 < /i < ^ in the 
remainder. We can write the outcome of the maximization 
problem in between brackets in ( [12] ) as 



r(p) = 




if < p < pth, 

( 1 + ;S^) if ^^- 

(13) 

Then we can rewrite the optimization problem in (TT2l in 
the same form as the point-to-point problem in (O with a rate 

3 The case with N± = N2 reduces to the single receiver problem. 
4 The time variable t is omitted for conciseness. 



function given in ( fl~3b . We next prove that this rate function 
is strictly concave. 

Lemma 4.1: The rate function r{p) in ( TT3b is a strictly 
concave function of power p. 

Proof: It is easy to show that r(p) is continuous, differ- 
entiable, and its derivative is decreasing with p; hence, it is a 
strictly concave function of p. ■ 

Now, based on this form of the optimization problem, we 
can directly use the results of Section [III] in the broadcast 
channel setting in order to identify the optimal transmission 
scheme for an energy harvesting transmitter. Note that as 
opposed to [ 1 3 1 and fl4l . our solution is valid for continuous 
energy arrivals as well as transmitters with various battery 
constraints. Once the optimal total transmit power over time 
is characterized, the power allocation among the users at each 
instant can be found using (13[ . 

V. Optimal Transmission Scheme with Battery 
Leakage 

In Sections Hill and [W] and references therein, the battery 
has been considered to be ideal, that is, there was no energy 
leakage. In this section, we consider the more realistic scenario 
of a battery that leaks part of the stored energy. 

The leakage rate of a battery depends on the type (Li-ion 
batteries have a lower leakage rate compared to the nickel- 
based ones), age and usage of the battery as well as the 
medium temperature. Moreover, even for a fixed type of 
battery and medium temperature, the leakage rate changes 
over time; the batteries leak most right after being charged. 
However, for simplicity, a constant rate leakage model is 
considered here. If the battery is non-empty at a given time 
instant, the energy is assumed to leak from the battery at a 
constant finite rate denoted by e > 0. Obviously no leakage 
occurs if the battery is empty. We use the same cumulative 
curve approach to model the battery leakage process. Note 
that the leakage rate e can alternatively be interpreted as the 
constant operation power of the node, that is, the circuit power 
needed to maintain the node awake. 

Definition 5.1 (Energy Leakage Curve): The energy leak- 
age curve L(t) is the amount of energy that has leaked from the 
battery in the time interval [0, t], t e M + , with L(0) = 0. Due 
to the constant leakage rate assumption, Lit) is a continuous, 
non-decreasing function whose right-derivative is given by 



L' + (t) = 



0, 



if E(t) < H(t) - L(t), 
otherwise. 



(14) 



To highlight the effect of leakage, we do not consider 
any minimum energy curve in this section, i.e., M(t) = 
Vt, and we focus only on discrete energy packet arrivals. 
Defining a maximum energy curve as U (t) = H(t) — L(t), the 
feasibility condition on the transmitted energy curve becomes 
< E(t) < U(t). We tackle again the problem of charac- 
terizing the feasible transmitted energy curve that transmits 
the most data over a given finite time interval [0, T]. The 
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Fig. 3. The Soo problem. 
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Fig. 4. The St problem. 

corresponding optimization problem can be stated as 

max T>(E(t)) = f n T r(E'(t))dt 
E(t)er 

such that < E(t) < U(t),Vt e [0,T]. 



(15) 
(16) 



Remark 5.1: Unlike the battery size constraint studied in 
Section [Hi] the battery leakage phenomenon does not translate 
into a minimum energy curve, but into a maximum energy 
curve obtained by removing the total leaked energy from the 
harvested energy curve. More importantly, the leakage curve is 
a function of the transmitted energy curve. Consequently, the 
maximum energy curve inherently depends on the transmitted 
energy curve, and hence, the solution framework presented in 
Section HII] does not directly extend to this setup. 

Throughout this section,we consider the discrete energy 
harvesting process in which the ?i-th energy packet of size 
E n arrives at time instant t n for n = 1,...,N. Without 
loss of generality, the first packet is assume to arrive at time 
t = (i.e., t\ — 0). We call this general setup the TV-packet 
problem. As before, we assume that the transmitter always has 
enough data to transmit. Below, we characterize the optimal 
transmission scheme first for the single-packet problem (i.e., 
N = 1), and then for the general TV -packet problem. 

A. The Single-Packet Problem 

We consider here the simplified problem consisting of a 
single energy packet E harvested at time t = 0. We refer to it 
as the single-packet problem. The solution of this problem will 
serve as a building block for the general TV-packet problem. 



First, let us treat the single-packet problem with infinit^f] 
deadline constraint (i.e., T = oo), and denote it by Soo- It 
is depicted in Fig. [3] Following Section [Till it is not hard 
to show that the optimal transmitted energy curve E(t) has 
to be piecewise linear, and the slope changes occur only if 
E(t) intersects U(t). Consequently, the optimal E(t) for the 
Soo problem is as shown in Fig. [3] the node transmits at a 
constant power p until the battery runs out of energy. One can 
see that there is a trade-off in the choice of p: while it is more 
energy efficient to transmit at lower power for a longer period 
of time, the longer the transmission time, the more energy 
will be wasted due to leakage. The optimization problem in 
([T5l>-([T6li becomes 



max T>(E(t)) = ^ r(p). 



p>0 



(17) 



Assuming that r{p) is a strictly concave increasing function 
with r(0) = 0, and a finite leakage rate e, the function 
f(p) — achieves its maximum at a finite p G R + , as 
shown in Appendix lAl We denote the corresponding optimal 
value by p*. Note that while the total amount of transmitted 
data is proportional to E, p* is independent of E. Summa- 
rizing, the optimal transmission strategy for the Soo problem 
is to transmit at constant power p* until the battery is empty. 
The total amount of transmitted data is 



r(p*). 



We next consider the single-packet problem with a fixed 
transmission deadline T, and denote it by St- It is de- 
picted in Fig. 2] and the following notations are defined: 
U(t) 4 H(t) - et (we assume U(t) > for all < t < T, 
as otherwise the problem is equivalent to the Soo problem). 
We denote the point (T, U(T)) by A. Finally, the slope of the 
line segment from the origin to A is denoted by s. We have 
s = E/T — e. As before p* denotes the value that maximizes 
the function f(p). Note that, as shown in Appendix lAl f(p) is 
strictly decreasing forp > p*. Hence, building on the solution 
derived for the Soo problem, the solution of the St is easily 
derived: 

> if s < p*, transmit at constant power p* until the battery 
is empty. 

> else, transmit at constant power s during the whole [0, T] 
interval. 

In short, the optimal transmission strategy for the St problem 
is to transmit at constant power p = max (p*,s) for a time 
duration -J^ (that is, until the battery is empty), and remain 
silent afterwards. The amount of transmitted data is ^p^r(p). 

B. The N -Packet Problem 

We consider here the general N packet problem with finite 
deadline constraint T, denoted as the Nt problem. 

We start with the following lemma which proves that the 
optimal solution of the Nt problem can be emulated in the 
equivalent St problem with E = ^n=i En - That ' s ' nav i n g 
all energy packets at time t = is at least as good as having 



5 Note that, in the case of energy leakage, potential transmit time is finite 
when the number of harvested energy packets is finite as the available energy 
decays to zero even if no data is transmitted. 



them arrive over time. Let us denote by D^ T and D$ T the 
optimal solutions (in terms of total transmitted data) of the 
Nt and equivalent St problems, respectively. 

Lemma 5.1: The optimal solution of the Nt problem can 
be obtained in the equivalent St problem with E = Yln=i ^n- 
That is, we have Ds T > Dn t . 

Proof: Consider the optimal curve for the Nt problem, 
and divide the [0, T] time interval into N sub-intervals: [ii, t 2 ], 
[t 2 ,t 3 ], [t n ,t n+1 ], [t N ,T]. We denote by T n the 
duration of the n th interval, i.e., T n = t n+ i — t n for 
n = 1, . . . , N — 1, and Tn = T — t^. From Theorem l3.3l we 
know that the optimal transmitted energy curve is a piecewise 
linear function, which is composed of constant power periods 
possibly separated by silent intervals (i.e. horizontal segments) 
in case the battery runs out of energy. Accordingly, we define 
the optimal solution of the Nt problem by the sequences 
{Pi,p 2 , ■ ■ ■ ,P N } and {T 1 ,T 2 ,.. L1 T N }, meaning that the 
node transmits for time T n (with T n < T n ) at power p n > 
in the n th interval. The node is silent in the remainder of the 
interval, i.e., for time T n — T n . 

The data transmitted by this transmission strategy 
is Djy T = E n =i r (Pn)- The total transmit energy is 
Y,n=i T n P n , while the total energy leakage is ^Y,n=i T n- 
Since the optimal solution should eventually empty the battery, 
we have 
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We now argue that this optimal solution can be emulated in 
the St problem with E = Y^,n=i -®n- Consider the following 
transmission strategy E(t) for the St problem: transmit at 
constant power equal to p 1 for time T\, followed by p 2 
for time T 2 , and so on, ending with p N for time Tn- 
By construction, this strategy transmits the same amount of 
data Dn t as the optimal solution of the Nt problem. We 
conclude the proof by showing that this strategy is feasible, 
that is, E(t) < U(t) for all t 6 [0,T]. Since the node is 
constantly transmitting during the interval [0, J2n=i Tn]> the 
curve U (t) is constantly decreasing during this interval at 
rate e, i.e. U{t) = J2n=i E n 
have 
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where the equality ( f20b follows from ( fT8l ). This proves the 
feasibility of E(t). 

Having proved the achievability of D^ T in the equivalent 
St problem, the inequality Ds T > Dn t naturally follows. ■ 

6 We assume U(t) > in the considered time interval, as otherwise the 
problem can be divided into equivalent subproblems. 
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(a) Original Nt problem (b) Original Nt problem 

Fig. 5. Counterexample for the equivalence of the Nt and St problems, with N = 2 energy packets. 



The counterpart of Lemma 15.11 in the other direction does 
not always hold, that is, the optimal solution of the equivalent 
St problem cannot always be emulated in the original Nt 
problem. A counterexample can indeed easily be constructed, 
which is depicted in Fig. Part (a) of the figure depicts a 
2-packet problem with Ex = E 2 = E, t\ = 0, and t 2 > T/2. 
In part (b) the equivalent St problem is depicted. Let the 
optimal transmission power for the St problem be given by 
p = 2yr — e. This solution cannot be emulated in the original 
Nt problem. In fact, as shown in part (a), the node cannot 
transmit a constant power p during the full [0, T] time interval 
as the battery runs out of energy at time T/2, and the node 
has to remain silent during the time interval [T/2, ti\. 

However, in the following lemma, we provide a sufficient 
condition for the counterpart of Lemma 15.11 to hold. For this, 
we define Ai as the point on the U(t) curve corresponding to 
the time instant t-t+i, for i = 1, 2, . . . , N — 1, and Am as the 
point corresponding to time t — T, as illustrated in Fig. [6] 

Lemma 5.2: If the line segment from the origin to the point 
Am does not cross the curve U(t) at any other point than 
{Ai, . . . iAm^, then the optimal solution of the St problem 
with E = 52n=i E n can b e obtained in the Nt problem, and 
Dn t > Ds T - This sufficient condition is expressed by the 
following N — 1 inequalities: 
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Proof: First note that the set of inequalities in ( l22l 
expresses that the line segments from origin to points 
Ax, Ai, . . . , Am-x have slope which are all greater than that 
of the segment from origin to Am- This requires that the line 
segments from the points A\, A 2 , . . . , Am-i to the point Am 
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Fig. 6. A Nt problem satisfying the conditions of Lemma 1531 



have slopes si,s 2 , ■ ■ ■ , sm-i respectively, which all are lower 
than or equal to the slope sq — n ^ — - 
from the origin to point Am- 



e of the segment 
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for i = 1, ...,N — 1. An illustration of an Nt problem 
satisfying the conditions of this lemma is given in Fig. [6] 

Consider now the St problem with E — z~2n=i E n- 
Remember that the optimal scheme for the St problem 
requires transmitting at constant power p for a duration of 
with p — max (p*, so) > sq. We now argue that this 
solution can be emulated in the Nt problem. Consider the 
following transmission strategy for the Nt problem: transmit 
at p whenever the battery is non-empty, and remain silent 
otherwise. By construction, this strategy is feasible. Again 
consider the N time intervals between energy arrivals [£1,^2], 



[t2j*3]> ■ • • > [tN, T] of durations T\, T2, . . . , Tn, respectively. 
We denote by T n (with T n < T„) the time for which the node 
is transmitting in the n th interval. The total transmission time 
is then given by T tot = XmLi Moreover, combining the 
inequalities in (l23l with the fact that p = max(p* , sq) > sq, 
we have that p > Sj for i = 1, . . . , N — 1. This ensures that 
the considered strategy uses up the whole available energy 
by time T, i.e., E(T) = U(T). Then, by the conservation of 
energy, the transmit and leakage energies must sum to the total 
harvested energy: 
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from which we get that T to t — -p^, just like for the optimal 
solution of the St problem. This transmission strategy thus 
transmits the same amount of data Ds T as the optimal solution 
of the St problem. Consequently, under the conditions given 
in the theorem, the inequality Dn t > Ds T holds. ■ 

Building on the two previous lemmas, the following theorem 
can be formulated. 

Theorem 5.3: If the inequalities in ( 1221 hold, then: 

(i) Dm t = T>s T , that is, the optimal solutions of the Nt 
problem and the St problem with E — Yl n =i are 
equivalent. 

(ii) The optimal transmission strategy for the Nt problem 
is to transmit at constant power p whenever the battery 
is non-empty, and remain silent otherwise, where the 
value p corresponds to the solution of the equivalent St 
problem: 




(25) 



The total amount of transmitted data is 



-(p). 

An illustration of the result in Theorem 15.31 is provided in 
Fig.|7]for N = 3 and p — p* > sq- Part (a) of the figure depicts 
the Nt problem, while its equivalent St problem is given in 
part (b). According to Theorem 15.31 for both problems the 
optimal strategy is to transmit at constant power p*. The only 
particularity of the Nt problem is the presence of silent zones 
in between energy packet arrivals. However, the distribution 
over time of these silent zones do not affect the total duration 
of transmission, guaranteeing the equivalence of both solutions 
in terms of amount of transmitted data. 

Now, building on Theorem 15.31 we can provide the optimal 
solution for any Nt problem. Consider all line segments 
connecting the origin to points Ai, i — 1,...,N. Among 
the segments that do not intersect U(t) other than at point 
{Ai, . . . , An}, we pick the one with the highest index, i.e., 
the rightmost end point. We denote this index by k. We 
can now consider the k first energy packets only, and solve 
the corresponding k packet problem with deadline Y^, n =i ^n, 
using the equivalence given in Theorem 15. 3 1 We then proceed 
recursively by considering the remaining N— k packet problem 
separately. This recursive algorithm is described next. It takes 



as inputs the number of packets N, the sizes of energy packets 
{En}n=i< anc l tne P a °ket interarrival times {T n }^ =1 . It returns 
as output the set of optimal transmission powers {p n }n—i' 
meaning that the optimal solution of the Nt problem is to 
transmit at constant power p n in the n th interval as long as 
the battery is non-empty. The optimality of the algorithm is 
proved in Appendix iBl 



Algorithm 5.1: NT-problem(N, E\,..., En, Ti, . . . , Tn) 
Input: 

• N: number of energy packets 

• {E n }n = i~. amount of energy in each packet 

• {T n }^ =1 : interarrival times 

Output: [p n }% =1 
Algorithm: 

1) Find the highest k G {1, . . . , N} such that 

Ei rp TT^k rp 

n=l ^ 2_^ n =\ ^ n /jQ) 

rri \ ^ k ryi 

n=l ln 2^in=l ln 

for all i S {1, . . . ,k — 1}. 



2) 



Pi = max P . ~ - e 



(27) 



for all i € {1, . . ., k}. 
3) If k < N, find the {Pn}n=k+i running 

N T -problem(N - k, E k+1 , E N , T k+1 , . . . , T N ) 



We conclude this section by identifying two special cases 
of the solution provided here: 

* The special case of an A^-packet problem without dead- 
line constraint can be solved by Algorithm |5T| by setting 
Tn = 00. In this case, the inequalities in (|26T > hold with 
k = N, and (|27| | reduces to pi — max(p*, 0) = p* for all 
i G {1, . . . , N}. Hence, the optimal transmission strategy 
for the A-packet problem without deadline constraint is 
to transmit at constant power p* whenever the battery is 
non-empty, and remain silent otherwise. 

> The special case of a perfect battery with no leakage 
is obtained by setting e = 0. In this case, p* — (as 
detailed in Appendix [A), and Algorithm 15.11 reduces to 
the solution proposed in ifTTI . 

VI. Conclusion 

We have considered a communication system with an 
energy harvesting transmitter. Taking into account various 
constraints on the battery we have optimized the transmission 
scheme in order to maximize the amount of data transmitted 
within a given transmission deadline. We have provided a 
general framework extending the previous work in [11 1 and 
[12 1 to the model with continuous energy arrival as well as 
time-varying battery size constraints. We have also showed that 
the proposed framework applies to the optimization of energy 
harvesting broadcast systems. Moreover we have studied the 
case of a battery suffering from energy leakage, for which 
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Fig. 7. Illustration of Theorem 15.31 with N = 3 energy packets and 



P = P 



the optimal transmission scheme has been characterized for a 
constant leakage rate. 

Appendix 

A. Properties of f(p) = 

Remember that r(p) is a non-negative strictly concave 
increasing function, with r(0) = 0. We prove here that the 
function f(p) = with p > 0, achieves its maximum at a 
finite p* 6 M + , and is strictly decreasing for p > p*. 

The derivative of f(p) is calculated as follows: 

r'(p)(p + e) -r(p) 



f'ip) 

We distinguish two cases: 



1) If e = 0, 



becomes 



f'ip) = 



(p + e) 2 



r'(p)p-r(p) 



(28) 



(29) 



which is analyzed as follows: 

• if p = 0, both the numerator and the denominator 
are zero. By l'Hopital's rule, we get lirxij,_>o f'(p) = 
r"(0)/2 < 0, which follows from the strict concavity 
of r(j>). 

• if p > 0, the numerator r'{p)p — r(p) is a strictly negative 
function. Indeed, the strict concavity of r(p) together with 
the fact that r(0) = guarantees that r(p) > r'(p)p for 
all p > 0. 

Overall, f'ip) is thus strictly negative for all p > 0. Hence, 
f(p) finds its maximum at p* = 0, and is strictly decreasing 
for p > 0. 



2) Consider now e > 0. The sign of ( I28i l is analyzed by 
focusing on its numerator only, which is rewritten for clarity 

as: 

nip) = r'(p)e + [r'{p)p — r(p)] (30) 
We analyze n(p) term by term: 



> The first term r'(p)e is a positive and strictly decreas- 
ing function, due to the increasing and strictly concave 
property of rip), respectively. 

> The term in between brackets r'(p)p — r{p) is equal to 
zero if p = 0, and a strictly negative (as shown above), 
strictly decreasing function for p > 0. Indeed, the strict 
concavity of r(p) guarantees that the derivative of this 
term r"(p)p is strictly negative for all p > 0. 

Consequently, overall nip) is a strictly decreasing function of 
p for p > 0. More precisely, the lower e, the more rapid 
the decrease of n(p) will be. The initial value at p = 
is positive and proportional to e: n(0) = r'(0)e > 0. On 
the other hand, the asymptotic value of n(p) as p — > oo is 
negative: linip^oo n(p) = linip^oo [r'(p)p - r(p)] < 0,where 
the inequality follows from the strict concavity of r(p) together 
with the fact that r(0) = 0. Between these two extremes, the 
strict decrease of nip) guarantees that f'(p) changes sign only 
once (from positive to negative) at some finite value denoted 
by p*, and that it will remain strictly negative for all p > p*. 
Hence, we have that: 

(i) f(p) has a unique maximum, which is achieved at some 
finite value of p > 0, denoted by p*. The lower the value 
of e, the lower the value of p* will be. 

(ii) f(p) is strictly decreasing for p > p* . 



B. Proof of Optimality of Algorithm 15.71 

If the inequalities in d26l) hold with k = N (as in Fig. |6]), 
the optimal solution provided in Theorem 15.31 is produced by 
Algorithm O in dZZ]) with k = N. 

Consider now that the inequalities in ( T26l i do not hold for 
k = N. Then, denote by k the highest k < N for which 
d26b holds. This situation is depicted in Fig. [8] We first argue 
that the optimal solution is such that it empties the battery 
before receiving the (k + l) th energy packet, i.e. before tk+i- 
Put differently, the optimal transmitted energy curve should 
intersect U(t) at a time t < tk+i- Assume that the opposite 
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Fig. 9. General Nt problem: suboptimal solution. 



holds, as depicted in Fig. [9] Then, at some time t' > tk+i, 
the slope of the transmitted energy curve E(t) has to increase 
in order to guarantee to empty the battery at time t = T 
(which is a necessary condition for optimality). However, it is 
easy to realize (see the dot-dashed curve in Fig. |9]l that such 
strategy is suboptimal since it violates Theorem 13. 31 Note that 
the feasibility of the dot-dashed curve in Fig. [9] in ensured by 
considering the largest k rather than any k satisfying d26b . 
Now, since the battery has to be emptied before receiving 
the (k + \) th energy packet, we can optimally decouple the 
problem. First, the k packet problem with deadline t^+i = 
2J n =i is solved independently. This subproblem satisfies 
the inequalities in d22b , such that Theorem 15.31 guarantees that 
its optimal solution is obtained by Algorithm |5.1| in d27l ). Then, 
proceeding recursively, the algorithm is run for the remaining 
N — k packet problem which can be considered as a new 
problem with an empty battery at the origin. 
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